Lattice study on nK scattering with moving wall source 



(N 

O 

(N 

a 

oo 



Ziwcn Fu 

Key Laboratory of Radiation Physics and Technology (Sichuan University), Ministry of Education; 
Institute of Nuclear Science and Technology, Sichuan University, Chengdu 610064, P- R- China. 

The s-wave pion-kaon ijiK) scattering lengths at zero momentum are calculated in lattice QCD 
with sufficiently light u/d quarks and strange quark at its physical value by the finite size formula. 
The light quark masses correspond to m-r: = 0.330 — 0.466 GeV. In the "Asqtad" improved staggered 
fermion formulation, we measure the nK four-point correlators for both isospin 1 = 1/2 and 3/2 
channels, and analyze the lattice simulation data at the next-to-leading order in the continuum three- 
flavor chiral perturbation theory, which enables us a simultaneous extrapolation of nK scattering 
lengths at physical point. We adopt a technique with the moving wall sources without gauge fixing 
to obtain the substantiable accuracy, moreover, for 1 = 1/2 channel, we employ the variational 
method to isolate the contamination from the excited states. Extrapolating to the physical point 
yields the scattering lengths as 771^03/2 = —0.0505(19) and = 0.1827(37) for / = 3/2 and 

1/2 channels, respectively. Our simulation results for -kK scattering lengths are in agreement with 
the experimental reports and theoretical predictions, and can be comparable with other lattice 
simulations. These simulations are carried out with MILC Nf = 2 -|- 1 fiavor gauge configurations 
at lattice spacing a ~ 0.15 fm. 

PACS numbers: 12.38.Gc 
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I. INTRODUCTION 

Pion-kaon (ttA') scattering at low energies is the sim- 
plest reactions including a strange quark, and it allows 
for an explicit exploration of the three-flavor structure 
of the low-energy hadronic interactions, which is not di- 
rectly probed in the tttt scattering. The measurement 
of nK scattering lengths is one of the cleanest processes 
and a decisive test for our understanding of the chiral 
SU(3) symmetry breaking of the quantum chromody- 
namics (QCD). In the present study, we will concentrate 
on the s-wave scattering lengths of "kK system, which 
have two isospin eigenchannels (/ = 3/2,1/2) in the 
isospin limit, and the low-energy interaction is repulsive 
for / = 3/2 channel, and attractive for / = 1/2 case, 
respectively. 

Experimentally, ttK scattering lengths arc obtained 
through ttK scattering phases using the Roy-Stainer 
equations. The experiments at low energies are an im- 
portant method in the study of the interactions among 
mesons [H-Ql: and these experiments have reported that 
the s-wave scattering length (oq) in the / = 3/2 chan- 
nel, TO7ra3/2 has a small negative value, namely, —0.13 ^ 
—0.05. Moreover, the on-going experiments proposed by 
the DIRAC collaboration Q to examine ttK atoms will 
provide the direct measurements or constraints on ttK 
scattering lengths. 

At present, theory predicts i:K scattering lengths with 
a precision of about 10%, and it will be significantly im- 
proved in the near future. Through the scalar form fac- 
tors in semi-lcptonic pseudo-scalar-to-pseudo-scalar de- 
cays, Flynn et al. Q extracted the t:K scattering length 
in the 7 = 1/2 channel as = -1-0. 179(17) (14). 

Three-flavor Chiral Perturbation Theory (xPT) has 
been used to predict the scattering lengths in the study of 
the low-energy i:K scattering, and small negative value 



was claimed as = —0.129 ^ —0.05. However, if 

the scattering hadrons contain strange quarks, xPT pre- 
dictions usually suffer from considerable corrections due 
to the chiral SU(3) flavor symmetry breaking, as com- 
pared with the case of the tttt scattering. Therefore, a 
lattice QCD calculation is needed to offer an alternative 
important consistent check of the validity of xPT in the 
presence of the strange quarks. 

To date, four lattice studies of tiK scattering length 
have been reported [l0l - [l3j . The first lattice calcula- 
tion of ttK scattering len gth in / = 3/2 channel was 
explored by Miao et al. [l^l using the quenched ap- 
proximation, and the value of was found to 
be —0.048. The first fully-dynamical calculation using 
Nf = 2+1 flavors of the Asqtad-improved [13, EBl stag- 
gered sea quark [l^, [l^l was carried out ll| to calculate 
the 7 = 3/2 scattering length for = 0.29 - 0.60 GeV, 
and further indirectly evaluate the 7 = 1/2 scattering 
length on the basis of xPT. They obtained a small nega- 
tive value of ~ —0.0574 for 7 = 3/2 channel and 
a positive value of TO7rQ'i/2 = 0.1725 for 7=1/2 channel, 
respectively. Nagata et al. fulfilled first direct lattice cal- 
culation on 7 = 1/2 channel [l3| using the quenched ap- 
proximation. They investigated all quark diagrams con- 
tributing to both isospin eigenstates, and found that the 
scattering amplitudes can be expressed as the combina- 
tions of only three diagrams in the isospin limit. This 
work greatly inspires us to study ttK scattering. How- 
ever, they did not observe the repulsive interaction even 
for 7 = 3/2 channel at their simulation points, and their 
lattice calculations are relatively cheaper. Sasaki et al. 
observed the correct repulsive interaction for 7 = 3/2 
channel and attractive for 7=1/2 case, and they ob- 
tained the scattering lengths of 771,^03/2 = —0.0500(68) 
and = 0.154(28) for the 7 = 3/2 and 1/2 chan- 
nels, respectively [l^. Moreover, to isolate the contam- 
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ination from the excited states, they construct a 2 x 2 
matrix of the time correlation function and diagonahze 
it [l^l, this method will guide us to study ttK scatter- 
ing for 7=1/2 channel in a correct manner. It should 
be stressed that, to reduce the computational cost, they 
employed a technique with a fixed kaon sink operator 
for the calculation of ttK scattering length for 7 = 1/2 
channel and then an exponential factor is introduced to 
drop the unnecessary t-dependence appearing due to the 
fixed kaon sink time p^ . In this work, we will improve 
this technique by using a "moving" wall source without 
gauge fixing where the exponential factor is not needed 
any more. Thus, there is no satisfactory direct lattice 
calculation for 7 = 1/2 channel until now. 

In the present study, we will use the MILC gauge con- 
figurations generated in the pres ence of TV/ = 2-1-1 flavors 
of the Asq tad improved [l3.ll5l| staggered dynamical sea 
quarks [13, [l3 to study the s-wave ttK scattering lengths 
for both 7 = 1/2 and 3/2 channels. Inspired by the ex- 
ploratory study of TTTT scattering for 7 = channel in 
Ref. [isj . we will adopt almost same technique but with 
moving kaon wall source operator without gauge fixing 
for 7=1/2 and 3/2 channels to obtain the reliable accu- 
racy. We calculated all the three diagrams categorized 
in Ref. [l^, and observed a clear signal of attraction 
for 7 = 1/2 channel and that of repulsion for 7 = 3/2 
case. Moreover, for 7 = 1/2 channel, we employ the 
variational method to isolate the contamination from the 
excited states. Most of all, we only used the lattice sim- 
ulation data of our measured ttK scattering lengths for 
both isospin eigenstates to simultaneously extrapolate to- 
ward the physical point using the continuum threc-fiavor 
xPT at the next-to-leading order. Our lattice simulation 
results of the scattering lengths for both isospin eigen- 
channcls are in accordance with the experimental reports 
and theoretical predictions, and can be comparable with 
other lattice simulations. 

This article is organized as follows. In Sec. |IT]we de- 
scribe the formalism for the calculation of irK scattering 
lengths including the Liischer's formula [l9l42l] | and our 
computational technique of the modified wall sources for 
the measurement of ttK four-point functions. In Sec. IIIII 
we will show the simulation parameters and our concrete 
lattice calculations. We will present our lattice simula- 
tion results in Sec. IIVI and arrive at our conclusions and 
outlooks in Sec. |Vl 



A. ttK four-point functions 

Let us consider the ttK scattering of one Nambu- 
Goldstone pion and one Nambu-Goldstone kaon in the 
Asqtad-improved staggered dynamical fermion formal- 
ism. Using operators 07r(a;i), 07r(a;3) for pions at points 
xi, X3, and operators Ok{x2)^ Oxixi) for kaons at points 
X2,X4^, respectively, with the pion and kaon interpolating 
field operators defined by 

Ojr+(x, i) = -d{x,t)-f^u{x,t) , 
O^o (x, = — [?l(x, <)75w(x, t) - d(x, t)75d(x, t)] , 
Oiio{^,t) = s(x, i)75cJ(x, , 

OkA^,^) = s(x,i)75u(x,<), (1) 
we then represent the ttK four-point functions as 

CTrK{x4,X3,X2,Xi) = {Ok{x4)0.^{x3)01^{x2)01{xi)) , 

(2) 

where (• ■ ■ ) represents the expectation value of the path 
integral, which we evaluate using the lattice QCD simu- 
lations. 

After summing over spatial coordinates xi , X2 , X3 and 
X4, we obtain the ttK four-point function in the zero- 
momentum state, 

C*7rx(i4,i3,i2,ii) = ^^^^C'7rK(a;4,a;3,a;2,xi) , 

XI Xo X3 X4 

(3) 

where xi = (xi,ii), X2 = (x2,i2), 2:3 = (x3,t3), and 
Xi = (x4, ti), and t stands for the time difference, namely, 

To avoid the complicated Fierz rearrangement of the 
quark lines, we choose the creation operators at the time 
slices which are different by one lattice time spacing as 
is suggested in Ref. namely, we select ti = 0,t2 = 
l,t3 = t, and = t + 1. In ttK system, there are two 
isospin eigenstates, namely, 7 = 3/2 and 7 = 1/2, we con- 
struct the irK operators for these isospin eigenchannels 
as 

0l^(t) = TT+m+it+l), (4) 
where 



II. METHOD OF MEASUREMENT 

In this section, we will briefly review the formulas of 
the s-wave scattering length from two-particle energy in 
a finite box, with emphasis on the formulae for isospin 
7=1/2 ttK system. Also we will present the detailed 
procedure for extracting the the energies of ttK system. 
Here we follow the original derivations and notations in 

Rcfs. [il,[ii,[2ii|. 



OKo{t) = ^OKo(x,t), OK+(t)=^OK+(x,i) 

X X 
X X 

If we assume that the u and d quarks have the same 
mass, only three diagrams contribute to ttK scattering 
amplitudes [l2j . The quark line diagrams contributing to 
ttK four-point function denoted in Eq. (jS]) are displayed 
in Figure [U labeling them as direct (D), crossed (C) and 
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(a) D (b) C (c) R 

FIG. 1: Diagrams contributing to iiK four-point functions. 
Short bars stand for wall sources. Open circles are sinks for 
local pion or kaon operators. The thinker lines represent the 
strange quark lines. 

rectangular (R), respectively. The direct and crossed di- 
agrams can be easily evaluated by constructing the corre- 
sponding four-point amplitudes for arbitrary values of the 
time slices and ^4 using only two wall sources placed at 
the fixed time slices ti and ^2- However, the rectangular 
diagram (R) requires another quark propagator connect- 
ing the time slices ^3 and t^, which make the reliable 
evaluation of this diagram extremely difficult. 



Sasaki et al. solve this problem through the technique 
with a fixed kaon sink operator to reduce the computa- 
tional cost [ij . Encouraged by the exploratory works of 
the TTTT scattering at / = channel in Refs. [H, H^l, simi- 
larly, we handle this problem by evaluating T quark prop- 
agators on a L'^ X T lattice, each propagator, which corre- 
sponds to a wall source at the time slice t = 0, - ■ ■ ,T— 1, 
is denoted by 

^A^^„"G^(n") = -^"'^ (-.*)' ^<i<T-l, (6) 

n" X 

where D is the quark matrix for the staggered Kogut- 
Susskind quark action. The combination of Gt{n) which 
we apply for t:K four-point functions is schematically 
shown in Figure [Tl where short bars stand for the posi- 
tion of wall source, open circles are sinks for local pion 
or kaon operators, and the thinker lines represent the 
strange quark lines. Likely, the subscript t in the quark 
propagator G represents the position of the wall source. 
D, C, and R, are schematically displayed in Figure [I] 
and we can also expressed them in terms of the quark 
propagators G, namely. 



C^(t4, ^3, t2,ti) 

(7*^(^4, ^3,^2,^1) 

G^{t4^, ^3, t2,ti) 



J2 J2 {^'^ Tr [gI (x3, t3)Gt, (X3, t3)Gt (x4, U)Gt, (x4, U) 

X3 X4 

(l^c Tr [gI (x3, t3)Gt, (X3, t3)Gl^ (x4, h)Gt, (x4, ^4) 

X3 X4 

Y (l^c Tr \gI {^2,h)Gu {^2,t2)Gl (x3, h)Gt, {^zM) 



(7) 



where daggers mean the conjugation by the even-odd par- 
ity (—1)" for the staggered Kogut-Susskind quark action, 
and Tr stands for the trace over the color index. The 
hermiticity properties of the propagator G are used to 
eliminate the factors of 7^ . 

For TiK rectangular diagram in Figure [Ijc), it cre- 
ates the gauge- variant noise [3 113 ■ O^.^ can reduce 
the noise by fixing gauge configurations to some gauge 
( e.g.. Coulomb gauge), and select a special wall source 
to emit only the Nambu-Goldstone pion [25[, however, 
the gauge non-invariant states may contaminate the ttK 
four-point function. Alternatively, we perform the gauge 
field average without gauge fixing since the gauge depen- 
dent fluctuations should neatly cancel out among the lat- 
tice configurations. Besides these cancelations, the sum- 
mation of the gauge- variant terms over the spatial sites 
of the wall source further suppresses the gauge- variant 
noise. In our current lattice simulation we found that 
this method works pretty well. 

All three diagrams in Figure [T] are needed to be cal- 
culated to study the ttK scattering in both 1 — 1/2 and 
1 = 3/2 channels. Three types of the propagators can be 



combined to construct the physical correlation functions 
for Trif states with definite isospin. As it is investigated 
in Ref. [l^, in the isospin limit, the -kK correlation func- 
tion for / = 3/2 and 1/2 channels can be expressed as 
the combinations of three diagrams, namely, 

cllHt)^{0llht)\0[lH^)) + ^NjC ^NfR, 
ciKht) ^ {oilht)\oilHo)) =D- NfG, (8) 

where the operator denoted in Eq. (|4]) creates a ttK 
state with total isospin / and the staggered-flavor factor 
Nf is inserted to correct for the flavor degrees of freedom 
of the Kogut-Susskind staggered fermion For the 

pion and kaon operators it is most natural to choose the 
one in the Nambu-Goldstone channel. This is the choice 
for our current study. 

To calculate the scattering lengths for hadron-hadron 
scattering on the lattice, or the scattering phase shifts in 
general, one usually resorts to Liischer's formula which 
relates the exact energy level of two hadron states in a 
finite box to the scattering phase shift in the continuum. 
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In the case of ttK scattering, the s-wave -kK scattering 
length in the continuum is defined by 



ao 



lim 



tanJo(fc) 



(9) 



k is the magnitude of the center-of-mass scattering mo- 
mentum which is related to the total energy of the 
'kK system with isospin / in a finite box of size L by 



EiK = Vr 



-HP. 



(10) 



where the m-„^ is the pion mass, and mx is the kaon mass. 
We can rewrite Eq. ([TU]) to an elegant form as 



E. 



K 



(11) 



In the absence of the interactions between the tt and K 
particles, A;/tan(So(fc) — > oo, and the energy levels occur 
at momenta k = 27rn/L, [n is a integer), corresponding 
to the single-particle modes in a cubic box. (5o(fc) is the 
s-wave scattering phase shift, which can be evaluated by 



the Liischer's finite size formula |19l. |21|. 



tan(5o(A;)\ ^ _ ( k^ 

^ ■ -Zoo I 1, ■ 



k J ttL V" (27r/L)^ 

where the zeta function Zoo(l; q^) is denoted by 



2oo(i;9') 



1 



1 



— ^ — cr 



(12) 



(13) 



here q = kLji^ii) is no longer an integer, and Zoo(l;9^) 
can be efficiently calculated by the method described 
in Ref. [2^. We also discussed this technique in Ap- 
pendix 1^1 where we extend this discussion to the case 
with the negative . In the case of the attractive inter- 
action, k? on the bound state has a negative value, there- 
fore fc is pure imaginary, and (5o(fc) is no longer physical 
scattering phase shift [l3- -Zoo(l, 9^), however, still have 
a real value even for this case, hence tan b^iki) jk obtained 
by Eq. is also real. If \k^\ is enough small, we can 
consider tan Jn(fc)/fc as the physical scattering length at 
Trfc threshold flSj . 

The energy E^^ of system with isospin / can be 
obtained from -kK four-point function denoted in Eq. 
with the large t. At large t these correlators will behave 

as IHlii 



CIk{^) = Z^K cosh 
(-l)*Z;^cosh 







T 
2" 




('- 


T 



(14) 



where E^^ is the energy of the lightest ttK state with 
isospin /. The terms alternating in sign are a pecu- 
liarity of the Kogut-Susskind formulation of the lattice 



fermions and correspond to the contributions from inter- 
mediate states with opposite parity [13, [l^ . The ellipsis 
suggests the contributions from excited states which are 
suppressed exponentially. 

We should bear in mind that, for the staggered Kogut- 
Susskind quark action, there are further complications 
in itself stemming from the non-degeneracy of pions and 
kaons in the Goldstone and other channels at a finite lat- 
tice spacing. Briefly speaking, the contributions of non- 
Nambu- Goldstone pions and kaons in the intermediate 
states is exponentially suppressed for large times due to 
their heavier masses compared to these of the Nambu- 
Goldstone pion and kaon [ll, [2^ Thus, we sup- 

pose that 7tK interpolator does not couple significantly 
to other ttK tastes, and neglect this systematic errors. 

In our concrete calculation, we calculated the pion 
mass and kaon mass ttik through the methods dis- 
cussed by the MILC collaboration in Refs. [1^ [s^l in our 
previous study [3l| . In this work we evaluate total energy 
EIx of ttK system with isospin / from Eq. ([Tl)) . 

In the current study we also evaluate the energy shift 
SEj ~ E^j^ — (m^ + itik) from the ratios 



a(o,i)Cx(i,t + i)' 



X = D,G, andi?. 



(15) 

where C7r(0, t) and C/f (1, t-|-l) are the pion and kaon two- 
point functions, respectively. Considering Eq. ([8]), we can 
write the amplitudes which project out the 7=1/2 and 
3/2 isospin eigenstates as 

Rj^.{t) = i?^(t) + \NfR(^{t) - ^NfR^^it) , 
Ri^.it) = R''{t)-NfR^{t). (16) 

Following the discussions in Ref. [l^l , we now then can 
extract the energy shift SEj from the ratios 



Ri{t) = Zie 



-SEit 



(17) 



where Zi stands for wave function factor, which is the 
ratio of two amplitudes from the ■nK four-point function 
and the square of the pion two-point correlator and the 
kaon two-point correlator, and the ellipsis indicates the 
terms suppressed exponentially. In Ri{t), some of the 
fluctuations which contribute to both the two-point and 
four-point correlation functions neatly cancel out, hence, 
improving the quality of the extraction of the energy shift 
as compared with what we can obtain from an analysis 
through the individual correlation functions [llj . 

For 7 = 3/2 channel, we can use Eq. ((Ti)) or Eq. (|17p to 
extract the energy shifts 5E. We have numerically com- 
pared the fltting values from two methods, and found 
well agreement within statistical errors. In fact, using 
Eq. PT)) to extract the energy shift 5E has been exten- 
sively employed for the study of ttK scattering at 7 = 3/2 
case in Ref. [11|. Hence, in this work, we will only present 
the energy shifts 5E calculated by Eq. (fTTl) . and then its 
corresponding scattering lengths. 
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On the other hand, for 7=1/2 channel, the presence 
of the kappa resonance is clearly shown in the low en- 
ergy [l^l, and therefore it should be necessary for us to 
separate the ground state contribution from the contam- 
ination stemming from the excited states to achieve the 
reliable scattering length as it is investigated in Ref. [l^l . 
For this purpose, we will construct a 2 x 2 correlation 
matrix of the time correlation function and diagonalize 
it to extract the energy of the ground state. 



B. Correlation matrix 

For 7=1/2 channel, to separate the contamination 
from the excited states, we construct a matrix of the 
time correlation function. 



C{t) 



{o\ol{t)OMom (o|ot(t)o.(o)|o) 



(18) 

where O,^ (t) is an interpolating operator for the k meson 
with zero momentum, and O^^Kit) is an interpolating 
operator for t:K system which is extensively discussed 
in section FlI Al The interpolating operator employed 
here is exactly the same as these in our previous stud- 
ies in Refs. |3ll - l33| . the notations adopted here are also 
the same, but to make this paper self-contained, all the 
necessary definitions will be also presented below. 



K sector 



In our previous studies [3lH33j, we have presented 
a detailed procedure to measure kappa correlator 
{0\k^ {t)K{0)\0) . To simulate the correct number of quark 
species, we use the fourth- root trick [sj], which auto- 
matically performs the transition from four tastes to one 
taste per flavor for staggered fermion at all orders. We 
employ an interpolation operator with isospin 7=1/2 
and = 0+ at the source and sink. 



Oix) 



1 



(19) 



where g is the indices of the taste replica, rir is the num- 
ber of the taste replicas, a is the color indices, and we 
omit the Dirac-Spinor index. The time slice correlator 
for the K meson in the zero momentum state can be eval- 
uated by 



;r E E (s'A^M'i^^t) u^(o,o).^(o,o)) 



x,a,b g,g' 



where 0,x are the spatial points of the k state at source 
and sink, respectively. After performing Wick contrac- 
tions of fermion fields, and summing over the taste index. 



for the light u quark Dirac operator M„ and the s quark 
Dirac operator M^, we obtain [3l| 

=E(-)"(Tr[AC'(x,t;0,0)A77i'(x,t;0,0)]) , 



(20) 

where Tr is the trace over the color index, and x = (x, t) 
is the lattice position. 

For the staggered quarks, the meson propagators have 
the generic single-particle form, 

Cit) = ^^e-™'* + A[{-lYe-<' + {t^Nt~t), 

i i 

(21) 

where the oscillating terms correspond to a particle with 
opposite parity. For k meson correlator, we consider only 
one mass with each parity in Eq. (j2ip . namely, in our 
concrete calculation, our operator is the state with spin- 
taste assignment 7(8)7 and its oscillating term with 7075(8) 
7o75 [HI- Thus, the k correlator was fit to the following 
physical model, 

C«(t)=6„e-"»* + 6^^(-)*e-^-A* + (t-^iV,_t), (22) 

where ^^id are two overlap factors. In Figure [51 
we clearly note this oscillating term. 

We should bear in mind that, for the staggered Kogut- 
Susskind quark action, our k interpolating operator cou- 
ples to various tastes as we examined the scalar ao and 
a mesons in Refs. [ssl [36j . where wc investigated two- 
pseudoscalar intermediates states (namely, bubble con- 
tribution). In Ref. [3^ . we investigated the extracted k 
masses with and without bubble contribution for kappa 
correlators. We found that there are only about 2% ^ 5% 
differences in masses, although the bubble contributions 
are dominant in the k correlators at large time region. 
Thus, in the current study, we can reasonably assume 
that the k interpolator docs not couple remarkably to 
other tastes, and ignore this systematic errors for the k 
sector I33I. 



2. Off-diagonal sector 

The calculations of the off-diagonal elements in correla- 
tion matrix C{t) in Eq. (HH), namely, {0\Olj^{t)O^{0)\0) 
and (OlOj. (t)C'Tii-(O) are exactly the same as these in our 
previous study for non-zero momenta in Ref. [s^ , the no- 
tations adopted here are also the same, but to make this 
paper self-contained, all the necessary definitions will be 
also presented below. 

To avoid the complicated Fierz rearrangement of the 
quark lines, we choose the creation operators at the time 
slices which are different by one lattice time spacing as 
suggested in Ref. [2^, namely, we select ii = 0,^2 = 1, 
and ~ t for ttK — > k three-point function, and choose 
ti = 0, t2 ~ t, and is = t + 1 for k — > t:K three-point 
function. 
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FIG. 2: Diagrams contributing to nK k and k — >■ nK 
three-point functions. Short bars stand for wall sources. The 
thinker lines represent the strange quark lines, (a) Quark 
contractions of nK hi, where open circle is sink for local 
kappa operator, (b) Quark contractions of k — >■ nK, where 
open circle is sink for local pion operator. 

The quark line diagrams contributing to the k — > -kK 
and ttA' — > k three-point function are plotted in Fig- 



ure [2ja) and Figure [2ljb), respectively, where short bars 
stand for the position of wall source, open circles are 
sinks for local pion or kaon operators, and the thinker 
lines represent the strange quark lines. Likely, the sub- 
script t in the quark propagator G represents the position 
of the wall source. 

The ttK — >■ K three-point function can be easily evalu- 
ated by constructing the corresponding three-point am- 
plitudes for arbitrary values of the time slice ^3 using only 
two wall sources placed at the fixed time slices ti and t2 ■ 
However, the calculation of k — > t:K three-point function 
is almost the same difficult as that of the rectangular di- 
agram for ttK four-point correlator function, since it re- 
quires another quark propagator connecting time slices 
t2 and ^3. The k — > ttK and nK — > k three-point func- 
tions are schematically shown in Figure [H and we can 
also express them in terms of the quark propagators G, 
namely. 



C'^K^re(i3i ^2, ^i) - 
Extraction of energies 



^ (Re Tr[Gt, i^3,h)Gl (xa, ^3)0^ (xi ,h)]) , 

X3,Xi 
X2,Xi 

I 



(23) 



Through calculating the matrix of correlation function 
C{t) denoted in Eq. (|T8)) . we can separate the ground 
state from first excited state in a clean way. It is very 
important to map out "avoided level crossings" between 
the K resonance and its decay products (namely, tt and 
K) in a finite box volume, because the first excited state 
is potentially close to the ground state. This makes the 
extraction of the ground state energy unfeasible if we only 
utilize a single exponential fit ansatz. Since we can not 
predict a priori whether our energy eigenvalues are near 
to the resonance region or not, we find it always safe 
in practice to adopt the correlation matrix to analyze 
our lattice simulation data for isospin I ~ 1/2 channel. 
To extract the ground state, we follow the variational 
method [l^l and construct a ratio of correlation function 
matrices as 



M{t,tR) = C{t)C-\tR), 



(24) 



with some reference time slice (20[, which is assumed 
to be large enough such that the contributions to cor- 
relation matrix M(t,tfj) from the excited states can be 
neglected, and the lowest two eigenstates dominate the 
correlation function. The two lowest energy levels can 
be extracted by a proper fit to two eigenvalues A„(t, t^) 
{n = 1,2) of matrix Mit^tn). Because here we work 
on the staggered fermions, and we can easily prove that 



^n{t,tR) {n = 1,2) behaves as [3S 



>'n{t,tR) = ^„cosh 
(-l)*S„cosh 



T 

-En {t-- 



(25) 



for a large t, which mean <C ifl < t ^ T/2 to sup- 
press the excited states and the unwanted thermal con- 
tributions. This equation explicitly contains an oscillat- 
ing term. For the current study, we are only interested 
in eigenvalue Ai(t,ti?,), here non-degenerate eigenvalues 
Ai(i,ii?,) > ^2{t,tR) are assumed. In practice, we found 
that the oscillating term in Xi{t,tfi) is not appreciable 
for some tji, we can also adopt following simple fitting 
model [SJ, 



(26) 



Ai(i,ifl) = Acosh{-E{t-T/2)) , 



and the difference between the fitting models of Eq. (PS)) 
and Eq. ([25)1 is small. However, to make our extracted 
ground energy E for isospin 7 = 1/2 channel always re- 
liable, in this work, we will present the ground energy E 
calculated by Eq. (j25|) , and then its corresponding s-wave 
scattering lengths. 
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III. LATTICE CALCULATION 
A. Simulation parameters 

We used the MILC lattices with 2 + 1 dynamical flavors 
of the Asqtad-iniproved staggered dynamical fermions, 
the detailed description of the simulation parameters can 
be found in Refs. |ll,[i3|- One thing we must stress that 
the MILC configurations are generated using the stag- 
gered formulation of lattice fermions [37l - [39j with rooted 
staggered sea quarks |30| which are hypercubic-smeared 
(HYP-smeared) |4CH43| | . As it was shown in Refs. [44l |45| 
that HYP-smearing gauge links significantly improves the 
chiral symmetry. 

Wc analyzed ttK four-point functions on the 0.15 fm 
MILC lattice ensemble of 450 16^ x 48 gauge config- 
urations with bare quark masses aniud ~ 0.0097 and 
ams = 0.0484 and bare gauge coupling — 6.572, 

which has a physical volume approximately 2.5 fm. The 
inverse lattice spacing = 1. 358^13 GeV [li,[i3- The 
mass of the dynamical strange quark is near to its phys- 
ical value, and the masses of the u and d quarks are 
degenerate. To avoid the contamination from pions and 
kaons propagating backward in time, periodic boundary 
condition is applied to the three spatial directions while 
in the temporal direction, Dirichlet boundary condition 
is imposed, which reduce the original time extent of 48 
down to 24, moreover, it avoids the "fake effects" dis- 
cussed in Rcf. jT3| . 

B. Sources for isospin 1—1/2 channel 

To calculate the ttK correlation functions, we use the 
standard conjugate gradient method (CO) to obtain the 
required matrix element of the inverse fermion matrix. 
The calculation of the correlation function for the rect- 
angular diagrams naturally requires us to compute the 
propagators on all the time slices = 0, • • • , T — 1 of 
both source and sink, which requires the calculation of 
48 separate propagators in our lattice simulations. After 
averaging the correlator over all 48 possible values, the 
statistics are greatly improved since we can put the pion 
source and kaon source at all possible time slices, namely, 
the correlator Cii(t) is calculated through 

CnW = {{7:K)it)i7:K)H0)) 

ts 

The best-effort to generate the propagators on all time 
slices enables us to obtain the correlators with high preci- 
sion, which is vital to extract the desired scattering phase 
shifts reliably. 

For each time slice, six fermion matrix inversions are 
required corresponding to the possible 3 color choices for 
the pion source and kaon source, respectively, and each 



inversion takes about one thousand iterations during the 
CG calculation. Therefore, all together we carry out 288 
inversions on a full QCD configuration. As shown follow, 
this large number of inversions, performed on 450 con- 
figurations, provides the substantial statistics needed to 
resolve the real parts of the / = 1/2 and 3/2 amplitudes 
with reliable accuracy. 

In the calculation of the off-diagonal correlator, C21 (t), 
the quark line contractions results in a three-point dia- 
gram. Since in this three-point diagram the pion field 
and kaon field are located at the source time slice ts, 
ts + 1, respectively. We calculate the off-diagonal corre- 
lator C2i(t) through 

C2l(t) = (/^(t)(^if)t(0)) 

= + (28) 

where, again, we sum the correlator over all time slices tg 
and average it. As for the second off-diagonal correlator 
Ci2(i), the pion field and kaon field are placed at the 
sink time slices ts +t and t + ts + 1, respectively, which 
make the computation of Ci2{t) difficult. However, using 
the relation Ci2{t) = C2i(t), we can obtain the matrix 
element C12 for free. As it is studied in Ref. [i^, since 
the sink and source operators are identical for a large 
number of configurations, C{t) is a Hermitian matrix. 
The K — >■ ttK component agrees with ttK — >■ k within 
the error, but the statistical errors of the matrix element 
C12 should be larger than that of matrix element C21 
for a large t. Therefore, in the following analyses we 
substitute matrix element C12 by the complex conjugate 
of matrix element C21, which is not only to save about 
20% computation time, but also significantly to reduce 
statistical errors. 

For the k correlator, C22{t), we have measured the 
point-to-point correlators with high precision in our pre- 
vious work [sij . Therefore, we can exploit the available 
propagators to construct the K-correlator 

C22W = ^5](«t(t + t,)A.(t,)), (29) 

where, also, we sum the correlator over all time slices ts 
and average it. 

One thing we must stress that, in the calculation of 
the correlator {{TTK){TTKy) , we make our best-efforts to 
reliably measure the rectangular diagram, since the other 
two diagrams are relatively easy to evaluate. We found 
that the rectangular diagram plays a major role in this 
correlator. Therefore, we get it properly for the ttK sec- 
tor for isospin 7 = 1/2 channel. 

In this work, we also measure two-point correlators for 
pion and kaon, namely, 

G^t) = (0|7rt(t)7r(0)|0), 

GK{t) = {0\K\t)K{0)\0) , (30) 

where the G^{t) is correlation function for the pion with 
zero momentum, and the Gxit) is correlation function 
for the kaon with zero momentum. 
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IV. SIMULATION RESULTS 



U I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I. 



In our previous work |47| . we have measured the pion 
and kaon point-to-point correlators. Using these corre- 
lators, we can precisely extract the pion mass (ttitt) and 
kaon mass [rriK) WA - which are summarized in Table ID 
Using the same method discussed in Ref. and the 
MILC code for calculating the pion decay constants f^r, 
we precisely extract the pion decay constants [49| . 
which are in agreement with the previous MILC deter- 
minations at this same lattice ensemble in Ref. [l7[ • Here 
we used all the 631 lattice configurations of this ensemble. 
We also recapitulated these fitted values in Table ID 



TABLE L Summary of the pion masses, kaon masses and 
the pion decay constants. The third and fourth blocks show 
pion masses and kaon masses in lattice units, respectively, and 
Column five shows the pion decay constants in lattice units. 
The second block gives pion masses in GeV, where the errors 
are estimated from both the error on the lattice spacing a and 
the statistical errors in Column three. 





m^(GeV) 




arriK 




0.00970 


0.334(6) 


0.2459(2) 


0.3962(2) 


0.12136(29) 


0.01067 


0.350(6) 


0.2575(2) 


0.3996(2) 


0.12264(34) 


0.01261 


0.379(7) 


0.2789(2) 


0.4066(2) 


0.12425(27) 


0.01358 


0.392(7) 


0.2890(2) 


0.4101(2) 


0.12482(32) 


0.01455 


0.406(7) 


0.2987(2) 


0.4134(2) 


0.12600(26) 


0.01940 


0.466(8) 


0.3430(2) 


0.4300(2) 


0.12979(27) 



A. Diagrams D, C, and R 

The ttK four-point functions are calculated with same 
lattice configurations using six u valence quarks, namely, 
am^ = 0.0097, 0.01067, 0.01261, 0.01358, 0.01455 and 
0.0194, where is the light valence u quark mass. They 
all have the same strange sea quark mass arris = 0.0426, 
which is fixed at its physical value [l^ . 

In Figure [3] the individual ratios, which are defined in 
Eq. ([T5|) corresponding to the diagrams in Figure [U 
{X = D,C and R) are displayed as the functions of t 
for arrix = 0.0097. We can note that diagram D makes 
the biggest contribution, then diagram C, and diagram R 
makes the smallest contribution. The calculation of the 
amplitudes for the rectangular diagram stands for our 
principal work. Clear signals observed up to f = 20 for 
the rectangular amplitude demonstrate that the method 
of the moving wall source without gauge fixing used here 
is practically applicable. 

The values of the direct amplitude R^ are quite close 
to unity, indicating that the interaction in this channel 
is very weak. The crossed amplitude, on the other hand, 
increases linearly, which implies a repulsion in / — 3/2 
channel. After an initial increase up to t '--^ 4, the rect- 
angular amplitude exhibits a roughly linear decrease up 



R°(t)-0.8 

<!>0<J><I><!><l><l><j><l><!><l><i>4><i><I><t><l><><I><I><l>^<I^ 



0.30 
0.25 
0.21 
0.15 
0.10 
0.05 
-0.01 
-0.C5 
-0.10 
-0.15 
-0.20 

2 4 6 8 10 12 14 16 18 20 22 24 
t 

FIG. 3: (color online). Individual amplitude ratios {t) 
of the TiK four-point function calculated through the mov- 
ing wall source without gauge fixing as the functions of t for 
arrix = 0.0097. Direct diagram (magenta diamonds) shifted 
by 0.8, crossed diagram (red octagons) and rectangular dia- 
gram (blue squares). 



iot^ 20, which suggests an attractive force between the 
pion and kaon in / = 1/2 channel. Furthermore, the mag- 
nitude of the slope is similar to that of the crossed ampli- 
tude but with opposite sign. These features are what we 
eagerly expected from the theoretical predictions @, H^] . 
We can observe that the crossed and rectangular ampli- 
tudes have the same value at t = 0, and the close values 
for small t. Because our analytical expressions in Eq. ([7]) 
for the two amplitudes coincide at t = 0, they should be- 
have similarly until the asymptotic t:K state is reached. 

In Figure m we display the ratio Ri{t) projected onto 
the isospin 7=1/2 and 3/2 channels for arrix = 0.0097, 
which are denoted in Eq. ([TB)l . A decrease of the ratio of 
Ri=3/2(t) indicates a positive energy shift and hence a re- 
pulsive interaction for 1 — 3/2 channel, while an increase 
of Ri=i/2{t) suggests a negative energy shift and hence an 
attraction for 1—1/2 channel. A dip at t = 3 for/ =1/2 
channel can be clearly observed [2^ . The systematically 
oscillating behavior for / = 1 /2 channel in the large time 
region is also clearly observed, which is a typical char- 
acteristic of the Kogut-Susskind formulation of lattice 
fcrmions and corresponds to the contributions from the 
intermediate states of the opposite parity [23, , this 
also clearly indicates the existence of the contaminations 
from other states rather than the pion-kaon scattering 
state Therefore, to isolate the potential contamina- 
tions, we will use the variational method |20j to analyze 
the lattice simulation data. As for 7 = 3/2 channel, this 
oscillating characteristic is not appreciable, we will use 
the traditional method, namely, using Eq. p7|) to com- 
pute the energy shift SE and then calculate the corre- 
sponding scattering length. 
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FIG. 4: (color online). Ri{t) for ttK four-point function at 
zero momenta calculated by the moving wall source without 
gauge fixing as the functions of t for arrix = 0.0097. Solid line 
in J = 3/2 is exponential fits for 7 < f < 16. 



B. Fitting analyses for 1 = 3/2 channel 

According to our discussions in Section|lIl in this work, 
we will make use of Eq. ([T7|) to extract the energy shift SE 
for 1 = 3/2 channel. Then we insert these energy shifts 
into the Eqs. © and ([T^ to obtain the corresponding 
s-wavc scattering lengths. Therefore, properly extract- 
ing the energy shifts is a crucial step to our final results 
in this paper. A convincing way to analyze our lattice 
simulation data is with the "effective energy shift" plots, 
a variant of the effective mass plots, where the propa- 
gators were fit with varying minimum fitting distances 
Dinin, and with the maximum distance Dmax either at 
the midpoint of the lattice or where the fractional statis- 
tical errors exceeded about 20% for two successive time 
slices. For each valence quark tUx, the effective energy 
shift plots as a function of minimum fitting distance D,nin 
for 7 = 3/2 channel are shown in Figure [5j The central 
value and uncertainty of each parameter was determined 
by the jackknife procedure over the ensemble of gauge 
configurations. 

The energy shifts a6E of ttK system for / = 3/2 
channel are extracted from the "effective energy shift" 
plots, and the energy shifts were selected by looking for 
a combination of a "plateau" in the energy as a function 
of the minimum distance Dmin, and a good confidence 
level (namely, x^) for the fit. We found that the effec- 
tive energy shifts for / = 3/2 channel have only relative 
small errors within broad minimum time distance region 
5 < D„iin < 10 and are taken to be quite reliable. 

We utilize the exponential physical fitting model in 
Eq. P?|) to extract the desired energy shifts for / = 3/2 
channel. In Figure|l]wc display the ratio Ri{t) projected 
onto the I = 1/2 and 3/2 channels for arrix = 0.0097, 



where we can watch the fitted functional form as com- 
pared with the lattice simulation data for 7 = 3/2 chan- 
nel. For the other five light u valence quarks, we ob- 
tain the similar results, therefore we do not show these 
ratio Ri{t) plots here. The fitted values of the energy 
shifts, 6Ei in lattice units and wave function factor Zj 
for 1 = 3/2 channel are summarized in Table |lTl The 
wave function Z factors are pretty close to unity and the 
X^/dof is quite small for / = 3/2 channel, indicating the 
values of the extracted scattering lengths are substan- 
tially reliable. 

TABLE II: Summary of the lattice simulation results for the 
energy shifts in lattice units for 7 = 3/2 channel. The third 
block shows the energy shifts in the lattice unit. Column four 
shows the wave function factors Z, Column five shows the 
time range for the chosen fit, and Column six shows the num- 
ber of degrees of freedom (dof) for the fit. All errors are 
calculated from jackknife. 

Isospin arrix aSE Z Range x^/dof 

0.00970 0.00621(53) 0.9880(59) 7- 16 0.0536/8 
0.01067 0.00615(52) 0.9893(58) 7- 16 0.0395/8 
/=§ 0.01261 0.00602(50) 0.9914(56) 7- 16 0.0226/8 
0.01358 0.00595(49) 0.9923(55) 7- 16 0.0176/8 
0.01455 0.00589(48) 0.9930(54) 7- 16 0.0142/8 
0.01940 0.00561(45) 0.9958(50) 7- 16 0.0067/8 



Now we can insert these energy shifts in Table ITll into 
the Eqs. ^ and to obtain the corresponding s-wave 
scattering lengths. The center-of-mass scattering mo- 
mentum fc^ in GeV calculated by Eq. (|lip . from which 
we can easily estimate its statistcal errors. Once we ob- 
tain the values of fc^ , the s-wave scattering lengths oq in 
lattice units can be obtained through Eqs. ^ and (|12p . 
All of these values are summarized in Table IIIII Here we 
utilize pion masses and kaon masses given in Table ID The 
errors of the center-of-mass scattering momentum k and 
the s-wave scattering lengths come from the statistical 
errors of the energy shifts energies SE, pion mass TOtt 
and kaon mass ttlk- 

TABLE III: Summary of the lattice simulation results of the s- 
wave scattering lengths for 1 = 3/2 channel. The third block 
shows the center-of-mass scattering momentum in GeV, 
Column four shows the s-wave scattering lengths in lattice 
units, and Column five shows the pion mass times scattering 
lengths. 



Isospin anix 


fc^[GeV2] 


do 




0.00970 


0.00350(27) 


-0.558(55) 


-0.137(13) 


0.01067 


0.00357(28) 


-0.569(55) 


-0.146(14) 


/ = 1 0.01261 


0.00366(30) 


-0.582(55) 


-0.162(15) 


0.01358 


0.00374(30) 


-0.593(56) 


-0.171(16) 


0.01455 


0.00379(34) 


-0.600(60) 


-0.179(18) 


0.01940 


0.00396(37) 


-0.624(62) 


-0.214(21) 
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FIG. 5: (color online). The effective irK energy shift plots, aSE as the functions of the minimum fitting distance Dmin in the 
fit for 1 = 3/2 channel. The effective ttK energy shift plots for 7 = 3/2 channel have only relative small errors within a broad 
minimum distance rGgion 5 ^ Dmin < 10. 



C. Fitting analyses for 1 = 1/2 channel 

In Figure [HI we show the real parts of the diagonal 
components {ttK — )■ ttK and k — k) and the real part 
of the off-diagonal component ttK — !• k of the correla- 
tion function C{t) denoted in Eq. (fTSl) . Since C{t) is 
a Hermitian matrix, we will substitute the off-diagonal 
component k — > ttK by ttK — >■ k to reduce statistical 
errors in the following analyses. 

We calculate two eigenvalues An(t, ii?) (n — 1,2) for 
the matrix M{t,tfj) in Eq. (|24p with the reference time 
tfl = 7. In this work, we are only interested in the eigen- 
value Xi{t,tj^) ^. In Figure [7] we plot our lattice results 
for Ai (t,tij) for each valence quark in a logarithmic 
scale as the functions of t together with a correlated fit to 
the asymptotic form given in Eq. (j25p . From these fits we 
then extract the energies that will be used to determine 
the s-wave scattering lengths. 

To extract the energies reliably, we must take two ma- 
jor sources of the systematic errors into consideration. 
One arises from the excited states which affect the cor- 
relator in low time slice region. The other one stems 
from the thermal contributions which distort the corre- 
lator in high time slice region. By denoting a fitting range 



^ In our previous study [33 
behavior of X2(t, tji). 



, we have preliminarily examined the 



10' 
10' 
10' 
10^ 
10^ 
10^ 
10' 
10" -i 
10' 
10-^ 
10-^ 
10^^ T 
10"' -i 
10 



-l — 1 — I — 1 — I — I — I — 1 — I — 1 — I — 1 — I — 1 — I — 1 — I — 1 — I — 1 — I — I — l- 



o 


G[k— >k] 





G[tcK->k] 


□ 


G[jiK->iiK] 



Go H R 



B H 







10 12 14 16 18 20 22 24 



FIG. 6: (color online). Real parts of the diagonal components 
{ttK — >■ nK and k — > k) and the real part of the off-diagonal 
component nK — !> k of the time correlation function C{t). 
Occasional points with negative central values for the diagonal 
component k — ^ k and the off-diagonal component nK k 
are not plotted. 



[imin,imax] and varying the values of the tmin and tmax, 
we can control these systematic errors. In our concrete 
fitting, we take tmin to be 4- 1 and increase the refer- 
ence time slice to suppress the excited state contam- 
inations. Moreover, we select imax to be sufficiently far 
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away from the time slice t = T/2 to avoid the unwanted 
thermal contributions. The fitting parameters tji, t,uijj 
and tinax are tabulated in Table IIVI The corresponding 
fitting results in the reasonable values of x^/dof. The 
X^/dof together with the fit results for the energies for 
the ground state aE are also listed in Table HVl 



TABLE IV: Summary of the lattice simulation results for the 
fitted values of the energy eigenvalues for the ground state for 
J = 1/2 channel. The third block shows the energy for the 
ground state in lattice units. In the table we list the reference 
time tR, the lower and upper bound of the fitting range, tmin 
and imax, the number of degrees of freedom (dof) for the fit 
quality x^/dof. All errors are calculated from jackknife. 



Isospin 



aE 



tn 



xVdof 



0.00970 
0.01067 
0.01261 
0.01358 
0.01455 
0.01940 



0.6260(6) 
0.6412(7) 
0.6694(10) 
0.6828(10) 
0.6952(11) 
0.7561(12) 



14 
14 
14 
14 
14 
14 



2.37/3 
2.83/3 
3.63/3 
4.55/3 
5.57/3 
13.1/3 



Now we can insert these energies in Table IIVI into 
the Eqs. and ([T^ to obtain the scattering lengths. 
The center-of-mass scattering momentum fc^ in GeV cal- 
culated by Eq. (|TT|) and thence the corresponding s- 
wave scattering lengths in lattice units obtained through 
Eqs. @ and are summarized in Table [V] Here we 
use the pion masses and kaon masses given in Table ID 



The errors of the center-of-mass scattering momentum 
k and the s-wave scattering lengths are calculated from 
the statistic errors of the energy shifts energies SE, pion 
mass rriTr and kaon mass niK- 



TABLE V: Summary of the lattice simulation results of the s- 
wave scattering lengths for / = 1/2 channel. The third block 
shows the center-of-mass scattering momentum in GeV, 
Column four shows the scattering lengths in lattice units, and 
Column five shows the pion mass times scattering lengths. 



Isospin anix 


fc^[GeV^] 


ao 


m^ao 


0.00970 


-0.00884(27) 


2.32(29) 


0.543(34) 


0.01067 


-0.00904(29) 


2.37(28) 


0.588(38) 


J = i 0.01261 


-0.00956(42) 


2.34(29) 


0.690(65) 


0.01358 


-0.01010(47) 


2.69(30) 


0.778(78) 


0.01455 


-0.01075(52) 


2.56(37) 


0.887(94) 


0.01940 


-0.01343(68) 


2.94(55) 


1.498(181) 



To examine the effects of the contaminations from the 
excited states for / = 1/2 channel, we denote the ratios of 
Coo(i) and EY[C~-^{tji) C(f)]oo which is the lowest eigen- 
value of C~^{tii) C{t), to the ttK correlator [l^. 



Roit) 
Do{t) 



Coo(i) 



Coo(iii) ' 

EV [C-\tn)C{t)]^ 



(31) 



In the upper panel of Figure [51 we show Ro(t) (magenta 
diamonds) and -Do(0 (blue octagons) at TOtt — 0.33 GeV. 
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FIG. 8: (color online). Ro{t) (magenta diamonds) and Do{t) 
(blue octagons) at (a) anix — 0.0097 or m,r — 0.33 GeV and 
(b) am^ = 0.0194 or m,r ^ 0.466 GeV for / = 1/2 channel. 



We can note that the difference of two ratios is small. 
This suggests that the contamination from the excited 
states is negligible at this light quark mass. However, 
from the bottom panel of Figure HI we observe that the 
contamination from the excited states is not negligible 
at ~ 0.466 GeV, and the diagonalization obviously 
changes the characteristics of the ratio, since the ttK in- 
terpolative operator for / = 1/2 channel has a large over- 
lap with the excited states [Ij]. Therefore, we further 
confirmed that the separation of the contamination from 
the excited states is absolutely necessary for the heavy 
quark masses [l3| when wc study the nK scattering for 
7=1/2 channel. 

Using the fitting models discussed in Rcf. [2^, we ex- 
tract the pion and kaon masses [131 . And using the fit- 
ting model in Eq (j^ . we calculate the kappa mass [s^ . 



In Figure [9l we display mK,TOif,m^, and ttK threshold 
TT^TT+K in lattice units as the functions of the pion mass 
TOtt- Wc observe that, as the valence quark mass in- 
creases, irK threshold grows faster than k mass and, as a 
consequence, for the heavy quark masses, t:K threshold 
is close to the k mass. In Figure 13 we can clearly note 
that ttK threshold is very close to the k mass for light 
quark mass am^ = 0.0097. This can be in part explained 
why the separation of the contamination from the excited 
states is indispensable for the large quark masses [l^ . 
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FIG. 9: (color online). Characteristics of m-K, m_fc, m-^ and 
rriTr+K in lattice units as the functions of the pion mass. 



D. Chiral extrapolations and scattering length 

In the present study, we employ a reasonable small pion 
masses mir, namely, mj!- ~ 0.330 — 0.466 GeV, which are 
still considerably larger than its physical ones, therefore, 
we need to extrapolate the s-wave nK scattering lengths 
toward the physical point. For this purpose, we employ 
the formula predicted by SU(3) chiral perturbation the- 
ory to next-to-leading order (NLO) 0, 1, [lil, [H HSl . In 
SU(3) chiral perturbation theory at NLO, we provide the 

1=1 1=1. 

continuum SU(3) xPT forms of Oq ^ and an ^ , which 

can be directly constructed from Refs. [sl. [TTI,]l3j . 



7=2 

2 



















-1 



32m^mx 



2 

J n 

32m7fTO/<- 
P 



^'^'^h (A 

J TT 



i5(Ax) + 



1 



16^2/2 



1 



167r2/2 



X./(A, 



,mK) 



(32) 



where we plugged in the pion mass m^, the kaon mass rized TableUl Lc,(A^) and L,rK{Ay^) = 2Li + 2L2 + — 
rriK, and the pion decay constant /jr, which are summa- 
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2^4 — j^fi/2 + 2£fi + L8 are low energy constants defined in 
Ref. [51[ at the chiral symmetry breaking scale A^. We 
should bear in mind that the expressions in Eq. (|32|) are 
written in terms of full /tt, and not its chiral limit value. 



The XnK (^x^'^-^'^k) and X^^^ (A^, m^r, mif ) are the 
known functions at NLO which clearly depend upon the 
chiral scale A^ with chiral logarithm terms, namely, 



X^K (Ax:W^i"iif) = K^ln— + In^^ + K,, In— + — "T-ifTOTT + Ktanarctan I — — _^ V 



^ A2 ' A2 ' 9 ^ " ' \ mK + 2m^ V 2toj^ 



\ / 1 ™^ , / , ™|- , / 1 "^^ , 86 3 (2{mK~m^) j niK 

XttK (^X'^^TT'"^^) = K^ln— + K^ln-— + K,, In— + — m^TO^ + -Ktanarctan ' ' 



A2 ' ^ A2 ' " A2 ' 9 2 V + V 2mA' - 



. arctan — tt: — — , (33) 



TTlK — 2to^ V 2mK + JTlir 



with 



K 



5mi 



2(m2 



kk = 



Kta 



67m^m7r — Sm'^mK + 23rnj^rn^ 

9 (ml- - 777,2 ) 

2Am-Km%; — brnxml + 28mj^m'^ — 9m^ 



18(?77|. 



ml) 



9 rriK 
nrriKml - 16m2,r772 - 



- 

107771 



2(m2 



%{m\ - 7772 ) 



24771,^777^ — brrifciTi^ 



18(7772,-7772) 
SmKmTj ^2777 1. — niKiriTr — ml 



9 



rriK 



(34) 



In this work, we did not measure the 77 mass (777^), al- 
teratively, we utilize the Gell-Mann-Okubo mass-relation 
to determine 77-mass. To improve the xPT fit, in princi- 
ple, we can include all the lattice simulation data of the 
■kK scattering lengths for both 1 — 1/2 and 3/2 channels 
to perform the simultaneous fitting. However, the fit 
with the data of the scattering lengths for / = 1/2 chan- 
nel in 77i7r > 0.392 GeV significantly increases x^/dof, so 
we only use the scattering lengths for / = 1/2 channel 
in ttZtt < 0.392 GeV. The fitting results of ttK scatter- 

I—'df'2 1—1/2 

ing lengths, rriT^aQ and niT^aQ arc plotted by the 
dotted lines as the functions of 7772 in Figure [TUl The 
dotted lines arc the chiral extrapolation of the scattering 
lengths for both isospin eigenstates. The fit parameters 
L.„K, L^, and the scattering lengths rriT^aQ at the physical 
points (namely, = 0.140 GeV, 777^ = 0.494 GeV) [HI 
are also summarized in Table IVIi where the chiral scale 
Apj, is taken as the physical rj mass, namely, A^^ = 0.548 
GeV [13 as it is done in Ref. [i^. The cyan diamond 
points in Figure [TUl shows the values of the physical scat- 
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FIG. 10: m^-dependence of the -kK scattering lengths m,rao 
for both 1 = 1/2 and 1 = 3/2 channels. The dotted lines give 
the SU(3) xPT predictions at NLO. The cyan diamond point 
indicate its physical values. 



tering lengths. From Figure [TUl we note that our lattice 
simulation results for 7 = 3/2 scattering length agrees 
well with the one-loop formula, while scattering lengths 
for 7=1/2 channel have a large error, and is in reason- 
able agreement with the SU(3) xPT at NLO. 



TABLE VI: The fitted s-wave scattering lengths ryi^^ao at the 
physical point {m-^ = 0.140 GeV, 771^= 0.494 GeV). The 
chiral scale A;^ is taken as the physical 77 mass. 



xVdof 10^ ■ L^K 10' ■ L5 



7= 



7= 



0.287/7 1.48(5) 1.06(14) -0.0505(19) 0.1827(37) 

The fitted value of the L5 is reasonable consistent 
with the value evaluated by PACS-CS collaboration (ssj . 
and is smaller than the corresponding result evaluated 
by MILC coUaboration ^ and NPLQCD coUabora- 
tion ll|. The fitted value of Lj^k is also smaller than 
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the result evaluated by NPLQCD collaboration [ll|. The 
s-wave TiK scattering lengths for both 7 = 1/2 and 3/2 
channels are in agreement with the other lattice stud- 
ies [IMl. 



V. SUMMARY AND OUTLOOK 

In the present study, we carried out a direct lattice 
QCD calculation of the s-wave ttK scattering lengths 
for both isospin / = 1/2 and 3/2 channels, where the 
rectangular diagram plays a crucial role, for the MILC 
"medium" coarse (a = 0.15 fm) lattice ensemble in the 
presence of 2 + 1 flavors of the Asqtad improved stag- 
gered dynamical sea quarks, generated by the MILC 
collaboration. We employed almost same technique in 
Ref. [3l but with the moving wall sources without gauge 
fixing to obtain the reliable precision. We calculated 
all the three diagrams categorized in Ref. [l^l, and ob- 
served a clear signal of the attraction for /= 1/2 chan- 
nel and that of repulsion for / = 3/2 channel, respec- 
tively. Moreover, for / = 1/2 channel, we employed 
the variational method to isolate the contamination from 
the excited states. We further confirmed that the sepa- 
ration of the contamination is absolutely necessary for 
the heavy quark masses [isj when we study ttK scat- 
tering in / ~ 1/2 channel. Simultaneously extrapolat- 
ing our lattice simulation data of the s-wave scattering 
lengths for both isospin eigenstates to the physical point 
gives the scattering lengths m^a3/2 = —0.0505(19) and 
= 0.1827(37) for the / = 3/2 and 1/2 channels, 
respectively, which are in accordance with the current 
theoretical predictions to one-loop levels and the present 
experimental reports, and can be comparable with the 
other lattice studies [l0-[l3|- 

A clear signal can be seen for long time separation 
range in the rectangular diagram of the -kK scattering. 
Reducing the noise by performing the calculation on a 
larger volume or the smaller pion mass could further im- 
prove the signal to noise ratio for the rectangular dia- 
gram, and therefore obtain better results for the scat- 
tering length in / = 1/2 channel (T^]. Moreover, the 
behavior near the chiral limit is strongly affected by the 
chiral logarithm term, so giving an evaluation without 
the long chiral extrapolation is highly desirable to ensure 
the convergence of the chiral expansion [l3|. Further- 
more, tan(5o(fc)/fc in the low- momentum limit must be 
evaluated by the systematic studies with the different vol- 
umes and boundary conditions [isj . For these purposes, 
we are planning a scries of lattice simulations on MILC 
coarse, fine, and superfine lattice ensembles with con- 
centrating on the lightest accessible values of the quark 
masses, namely, in < 300 MeV. We anticipate that 
these future tasks should make the calculation of the rect- 
angular diagram more reliable. 

It is well-known that -kK scattering at / = 1/2 channel 
is a more challenging and interesting channel phenomeno- 
logically due to the existence of kappa resonance. The 
study of the s-wave -kK scattering at zero momentum is 



just first step in the study of the hadron interactions in- 
cluding s-quarks. However, it is particularly encouraging 
that ttJsT scattering for / = 1/2 channel can be reliably 
calculated by the moving wall sources without gauge fix- 
ing in spite of the essential difficulties of the calculation of 
the four-point functions, especially rectangular diagram. 
It raises a prospect that this technique can be successfully 
employed to investigate the k resonance, and so on. 

In our previous work [3l|, we have precisely evaluated 
the K mass, and found that the decay k — >■ ttK is only al- 
lowed kinematically for enough small u quark mass. This 
work and our preliminary lattice simulation reported here 
for TiK scattering lengths will encourage the researchers 
to study the k resonance. We are beginning a series of 
lattice investigations on the k resonance parameters with 
isospin representation of (/, Iz) — (1/2, 1/2), and the pre- 
liminary results are already reported in Refs. [33l |54|. 
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Appendix A: The calculation method of zeta 
function 

In this appendix we briefiy discuss one method for 
numerical evaluation of zeta function Zoo(s;(7^) in the 
center-of-mass system for any value of q^. Here we fol- 
low the original derivations and notations in Ref. f26j . 

The definition of zeta function Zoo(s; q^) in Eq. ([TH]) is 



47r-Zoo(s;g2)= ^ 



1 



(n2 



(Al) 



The zeta function 2oo(s;q2) takes a finite value for 
Res > 3/2, and Zoo(l;9^), is defined by the analytic 
continuation from the region Res > 3/2. 

First we consider the case of > 0, and we separate 
the summation in Zqq{s; q^) into two parts as 



nG 



[n? — q^Y 



(A2) 



The second term can be written in an integral form. 
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y ' = — y 

»1 



1 



ns) Jo 



1 



E 



(A3) 



The second term neatly cancels out the first term in 
Eq. ( Next we rewrite the first term in Eq. (|X3)) 

by the Poisson's resummation formula as, 



1 f' 



r{s) 



f dt f-^e"^' 

Jo .^^n,: 



^ ^'dtt-ie^-^^r^^'^'^e-^^-V*. (A4) 



r(s) Jo 



The divergence at s = 1 comes from the n = part of 
the integrand on the right-hand side, therefore we divide 



the integrand into a divergent part (n = 0) and a finite 
part (n ^ 0). The divergent part can be evaluated for 
Re s > 3/2 as 



1=0 ' 



The right hand side of this equation takes a finite value 
at s = 1. 

After gathering all terms we obtain the representation 
of the zeta function in the center-of-mass system at s = 1, 



J 



neZ3 ^ 1=0 ' neZ3 



where X^Lez^ stands for a summation without n = 0. 

For the case of < 0, it is not necessary for us to 
separate the summation in Zqq{s\ g^), and it can be also 
written in an integral form. 



y ' 



ne 



dt t'-^e""' 



E 



nez 



(in? — q^Y 



(A7) 



Following the same procedures, we arrive at the same 
expression in Eq. (jA6p . Hence, Eq. (|A6[) can be applied 
for both cases. 

We also note that, for negative g^, an exponentially 
convergent expression of the zeta function 2oo(s; q^) has 
been derived in Ref. (ssj . We numerically compared this 
representation of the zeta functions with that of above 
described representation, and found agreement. There- 
fore, we use Eq. ()A6P in this work. 
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